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Q-i, The weak AA" — > NN transition is studied in the valence quark model approach. The 

momentum transfer for this transition is so large that the short- distance two baryon dy- 
namics must be taken into account. The two baryon system is described in the quark 
cluster model and the weak transition amplitude is calculated by evaluating the matrix 
elements of the effective weak AS = 1 hamiltonian. The results indicate some qualitative 
differences when compared with those in conventional meson-exchange calculations. Es- 
pecially, we conclude that contributions of the AI = | transition are significant and that 
the discrepancy in the n — p ratio between theory and experiment could be resolved by 
including the direct-quark processes. 
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1. Introduction 



The hyperon A decays weakly into a nucleon and a pion in the free space. Two isospin 
modes, pn~ and nir°, share 64% and 36% of the total decay. If the decay goes through 
a AI = | vertex, the p-rr~ /mr ratio would be two to one except for a small correction 
due to the phase space difference. The experimental ratio is very close to the AI = | 
prediction, and thus support the AI = \ hypothesis for the hadronic weak decay. 

In the nuclear medium, the A — > A% decay is suppressed by the Pauli blocking on 
the final nucleon state, whose momentum is less than 100 MeV/c for the A decay at 
rest. Indeed, in heavy hypernuclei, the decay is predominantly the nonmesonic one, that 
is, AA" — > NN. If we assume that the initial A and the nucleon are at rest, then the 
final relative momentum of A^A" is about 420 MeV/c and thus is well above the Fermi 
momentum. 

The purpose of this report is to study the direct quark processes in the two-body 
AA" — > NN weak decays and to show qualitative differences from the conventional picture 
employing the meson exchange mechanism. We present a possibility to solve the problems 
that the meson exchange mechanism encounters. 

The nonmesonic decays of hypernuclei seem quite useful in studying the low energy 
nonleptonic weak interactions among quarks. The final relative momentum of AW is high 
enough to look into the short-distance component of the two nucleon system. This type of 
the hyperon decay may reveal a new aspect of the weak interaction under the influence of 
the strong interaction. An advantage of using the hypernuclear decay is that the process 
is selective in isospin, spin and orbital angular momentum for appropriate initial and final 
states of the hypernucleus. 

Theoretical study of the nonmesonic decay of hypernuclei has traditionally employed 
the meson (ir, K, p, etc.) exchange mechanism, where one of the meson-baryon vertices 
involves the weak transition s — > Accumulating experimental data, however, have 
revealed some difficulties in the meson-exchange picture. For instance, the so-called n — p 
ratio, i.e., the ratio R np of An nn v.s. Ap — > np decay in the nucleus, is predicted 
very small, R np ~ 0.1 — 0.4 in the meson-exchange picture. This is due to the strong 
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contribution of the tensor force, which is preferred at the large momentum transfer. The 
tensor force selects the S = 1, I = pn final state and therefore R np becomes small. The 
experimental data seem not to agree with the prediction, i.e., ~ 1 in decays of light 
hypernuclei. We argue that the direct quark process, which does not follow the 1 = 
selection rule, may enhance the n — p ratio. 

The mesonic weak decays of hyperons have been tested for the AI = | rule and are 
known to satisfy the rule to about 5% error. The same rule for the nonmesonic weak 
processes, like AiV — > NN, is not confirmed yet. Indeed, an analysis of the decay of the 
A = 3 and 4 hypernuclei claims that the AI = ~ rule may be violated 0. It is therefore 
urgent to clarify the mechanism of the AI = | rule in the free hyperon decays and to 
study whether the same mechanism restricts the nonmesonic decays to AI = | as well. 

In the study of the meson-exchange processes, the AI = ~ rule is assumed from the 
beginning, implemented in the A — > Ntt vertex. We instead employ the effective quark- 
quark weak hamiltonian, which contains both the AI = | and AI = | components. 
Although the A J = | part has a small overall coefficient, we will see that its matrix 
elements for the AA^ — > NN decay may not be small compared to the AI = | component. 

The paper is organized as follows. Sect 2 is devoted to the basic formulation of the 
present calculation. The effective hamiltonian and the quark cluster model wave functions 
are presented. Various approximations employed in this calculation are examined. We 
present the results of the calculation in sect. 3. They are used for studying qualitative 
differences between the direct-quark mechanism and the conventional meson-exchange 
picture. A conclusion is given in sect. 4. 

2. Formulation 

2.1. Effective Weak Hamiltonian 

The effective weak hamiltonian describing AS = 1 processes has been calculated by 
several author s||,f|,f|. It can be computed by analyzing the correction due to the strong 
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interaction on the pure weak vertex su — > du: 

G f 

if (purely weak: AS* = 1) = --J= sin6» c {u a s a ) V - A (dpUp) V - A (1) 

where 

{u a s a ) v ^ A = {u a ^{l - 75) s a ) etc. (2) 

and a and f3 denote the color of quarks and the color sum is always assumed. Note that 
there exists no strangeness-changing neutral current in the standard electro-weak theory, 
and thus the vertex is only for the left-handed quarks. 

This purely weak four-quark vertex clearly contains both the AI = | and AI = | 
components. It was pointed out, however, that the correction due to the strong interaction 
enhances the AI = ~ component while the AI = | is suppressed at the same time||. The 
strong correction is treated perturbatively at the scale Q 2 ~ M^, and only the lowest 
order diagrams are taken into account. The mechanism of the AI = | enhancement 
can be understood by realizing the anomalous dimensions of the AI = | and AI = | 
components in eq.(|l|) have opposite signs with each other. The AI = | (AI = |) 
anomalous dimension is positive (negative) and therefore, when the renormalization scale 
is moved down from the W boson mass to the low-energy hadronic scale (<1 GeV), 
the operator with the positive anomalous dimension is enhanced and vice versa. The 
renormalization group equation also induces new four-quark operators through operator 
mixings. Another contribution comes from the so-called penguin diagrams. They are 
purely AI = | and thus help the AI = | rule. 

Taking these effects into account, the low energy effective weak hamiltonian has been 
derived 0: 

HifiT 1 (Q 2 - V 2 ) = -% E K r O r (3) 
where the four-quark operators, Ok (k = 1, 2, 3, 5 and 6) are defined by[[| 

01 = (d a S a ) V - A (upUf3) V -A - (u a s a )v-A(dpUp)v-A (4) 

02 = {d a S a ) V - A {upU(3) V -A + (u a S a ) V - A (d(5Uf3)v-A 
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+ 2(d a s a )v-A(dpdp) V - A + 2(d a s a ) V - A {spSf3)v-A (5) 
3 = 3 {AI = \) + 3 {AI = \) (6) 

3 ( AI — - (d a S a ) V - A {upUf3)v-A + {u a S a ) V ~A{dpUf3)v-A 

+ 2{d a S a )v-A{df3dp)v~A — 3{daS a )v-A{spSp)v-A 

3 5 r 

3 ( A / = -) = - y(d a S a ) V - A (upUp) V-A + {u a Sa)v-A{df3Up)v-A 
— {d a S a )v-A{dpdi3)v-A 

5 = {d a S a ) V _ A (u^Up + dpdp + SpSp) V+A (7) 

6 = (d a sp) V -A{upu a + d p d a + spS a ) v+A (8) 

Among these operators, 3 contains a part that induces the AI = | transition, 3 (AI = 
|), while the others are purely AI = |. The coefficients can be calculated by solving 
the renormalization group equation to the one-loop QCD corrections. They depend on 
the mass of the top quark mt through the penguin diagrams, which generate the operators 
O5 and Oq with the V + A coupling. We find that the final results are insensitive to the 
choice of m t and here choose m t = 200 GeV/c 2 . The coefficients also depend on the energy 
scale fi 2 for the effective hamiltonian. In the present calculation, we choose two sets of 
the values given in ref. 0: fi = 0.24 GeV and 0.71 GeV. They are chosen so as to give 
a s (/i 2 ) = 1 for the QCD Aq C d parameter, Aq C d = 0.1 GeV and 0.316 GeV, respectively. 
The values of the coefficients K r used in the present calculation are given in Table 1. One 
sees that the two choices are not much different except for K§ and K$. We will find that 
the differences in the transition matrix elements for these choices are at most 10%. 

The most prominent feature of this effective hamiltonian is that the QCD correction 
enhances the 0\ component while the other terms are suppressed. This is the main 
mechanism for the AI = ~ enhancement as is explained above. Later we will compare 
the results with and without 3 (AI = |) in order to study the AI = | contribution. 

This effective hamiltonian has been used for the calculations of the nonleptonic decay 
of strange mesons and baryons[§]]. It is found that although the AI = | transition is indeed 
enhanced in those decays, the enhancement is not enough to account for the experimental 
data quantitatively. It was suggested^] that an additional AI = | enhancement arises 



Table 1: Two choices of strengths of the weak effective four-fermi vertices, taken from 
ref.[5]. We use the version with flavor-number dependent A and m t = 200GeV. The values 
of the CKM matrix elements are taken as the central values of those given in ref.[6]. 





fi (GeV) 


A< 4 > (GeV) 


K x 


K 2 


K 3 


K 5 


K 6 


I 


0.24 


0.10 


-0.284 


0.009 


0.026 


0.004 


-0.021 


II 


0.71 


0.316 


-0.270 


0.011 


0.027 


0.002 


-0.010 



from the mesonic correction in the chiral effective theory. It was also suggested that the 
decay amplitudes may be sensitive to the meson and baryon wave functions ||. 

2.2. Six-Quark Wave Function 

In calculating the decay amplitude for AA" — > NN, we employ the constituent quark 
model, which describes the spin-flavor structure of the ground-state baryons very well. 
Two-baryon systems are expressed by the quark-cluster-model wave functions [PUlPl- First, 
we assume that the baryon consists of three valence quarks, whose orbital wave function 
is a harmonic oscillator eigenstate, 

where £'s are the Jacobi coodinates and the Gaussian parameter b is chosen b = 0.5 fm. 
We here neglect the asymmetry due to the mass difference of the s and u quarks in the A 
wave function. The six quark wave functions are given by 

\AN) — A 6 \(p(l, 2, 3)0(4, 5, 6)xo(R)) 

\NN) = A 6 \<P(1, 2, 3)0(4, 5, 6) X (R)) (10) 

where A 6 is the antisymmetrization operator for six quarks, is the internal wave function 
of the baryon, and R is the relative coodinate of two baryons. Xo{R) (x(R)) is the initial 
(final) relative wave function. The flavor-spin part of is taken to be purely the SU(6) 
wave function. 
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In the present calculation, we choose the simplest relative wave functions, i.e., a 
Gaussian for the initial state and the plane wave for the final state. 




X (R)=exp{ik-R} (12) 



The relative momentum of the final state is determined by the realistic Q value AE = 
Ma — Mjv of the decay: k = 416 MeV/c. The Gaussian B parameter of the initial 
state is to be determined by the AN wave function in the hypernucleus. Suppose that 
we choose (unrealistically) B = yz/3 6, then the initial wave function is reduced to the 
(Os) 6 configuration in the harmonic oscillator shell model. This can be interpreted as a 
dibaryon state in which the initial A and iV are on top of each other (R — 0). Later we 
show the results of the calculation where we employ B = 1.838 fm (= v^x 1.3), that 
corresponds to a AiV system in the 4 He nucleus. 

Although these choices of the wave functions may not be totally realistic, they will 
clarify the qualitative difference between the meson-exchange and the direct-quark pro- 
cesses, which is the purpose of this study. An advanced calculation using the realistic 
two-baryon wave functions is under way. 

Because we employ the nonrelativistic valence quark picture for the wave functions, 
the effective hamiltonian is also approximated by adopting the Breit-Fermi nonrelativistic 
expansion up to p/m. Then the spin-orbital part of the operator contains terms which 
conserve parity, 

lxl, (<7j • Bj) x 1 

and those which break parity, 

(Bi ± Bj) ■ q i:j , (Bi ± Bj) ■ (Pi ± Pj) , % (Bi x Bj) ■ cfij, i x Bj) ■ (Pi ± P^ 

where qlj = — Pi and Pi = \{Pi + with pi (p^) the initial (final) momentum of the 
2-th quark. In this expansion, we take account of the SU(3) breaking effects due to the 
quark mass differences. Then the coefficients for these operators depend on the light 
quark constituent mass, m q = m u = 313 MeV and the mass ratio, m u /m s = 0.6. 
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Table 2: Possible initial and final quantum numbers for the initial L = transition 

channel isospin spin-orbital 
1 pA — > pn 

2 
3 
4 
5 
6 
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nA - 




%- 


-> 1 S'o 




8 










K 


9 






3 5x- 


- 3 P1 


fn 



So ~ 


-> So 


a v 




- 3 Po 


b p 








3 ^- 


- 3 ,Di 


dp 


3 ^- 




Bp 


3 5x- 


- 3 Pi 


fv 



In the present study, we restrict our initial state to L — 0. Table 2 shows nine possible 
combinations of L, S, J, and / for the initial and final states. We note that the 1 = 1 
final states are allowed both for (An — > nn) and (Ap — > pn), while the 1 = states are not 
possible for (An — > nn) . Thus we have 6 (Ap — > pn) and 3 (An — - > nn) matrix elements, 



which are labeled from a through / in Table 2, according to the widely used notationfll]. 
Among them, the channels 1, 3, 4, and 7 (a, c and d) are the parity conserving transitions, 
while the others violate the parity invariance. The amplitude d is not zero only for the 
tensor component of the weak interaction, which we neglect in our quark model calculation 
by truncating the p/m expansion. 

The transition amplitudes are calculated according to the standard quark cluster model 
approach. Remember that we take into account the full antisymmetrization among six 
quarks. One needs to calculate the exchange matrix elements as well as the direct ones. 

3. Results 

The results of the calculation are summarized in Table 3. Nine amplitudes give all 
the information for the AN — > NN weak decay from L = 0. We compare the results for 
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Table 3: Calculated transition matrix elements in 10 10 MeV l / 2 . The numbers without 
parenthesis (in parenthesis) are the results for the parameter I (II). 



channel 


full 


AI = 


| omitted 


OPE 


a p 


-6.66 


(-6.83) 


-0.02 


(0.04) 


4.52 


b p 


5.79 


(6.30) 


0.06 


(0.37) 


-24.4 


Cp 


2.70 


(2.48) 


2.70 


(2.48) 


4.52 


dp 





(0) 





(0) 


-83.0 


Cp 


-2.21 


(-2.01) 


-2.21 


(-2.01) 


-42.3 


fp 


-5.57 


(-5.02) 


-5.39 


(-4.83) 


19.9 


a n 


4.67 


(4.92) 


-0.03 


(0.06) 


6.39 


b n 


-3.96 


(-3.67) 


0.09 


(0.52) 


-34.5 


fn 


-7.49 


(-6.70) 


-7.62 


(-6.83) 


28.2 



the two choices of the weak hamiltonian parameters, which differ in the energy scale or 
Aqcd, given in Table 1 and see that their differences are small. The numbers given under 
the 11 AI = | omitted" are the results without the AI = | component of the 3 operator. 
In this case, the ratio a n /a p , b n /b p , and f n / f p are equal to \/2. The other amplitudes, 



c, d and e do not contain any AI = § component, because the final NN states have 
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1 = 0. We find that the amplitudes a and b get significant contributions from the AI = § 
component, while / has only a small contribution of AI = |. Thus we conclude that the 
AI = | transition can be studied in the weak decay starting from the AN 1 Sq state. 

In Table 4, we summarize the calculated decay rates with the initial spin averaged and 
the final states summed up. 

T P = ^\[al + bl + 3(cl + dl + el + a (13) 

r " = Wi [a " + 6 " + 3/ " ] (14) 

We again find that the AI = | component of O3 changes the total decay rate by as much 
as a factor two for the direct quark processes. 
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Table 4: Calculated observables. The numbers without parenthesis (in parenthesis) are 
the results for the parameter I (II). 





full 


AI = 


| omitted 


OPE 


T p (10 7 sec" 1 ) 


0.39 


(0.36) 


0.23 


(0.19) 


52.4 


T„ (10 7 sec" 1 ) 


0.39 


(0.32) 


0.33 


(0.26) 


6.8 


Rnp 


0.99 


(0.89) 


1.41 


(1.39) 


0.13 


V P 


2.13 


(1.96) 


4.64 


(4.47) 


0.34 


Vn 


8.43 


(6.12) 


2.23-10 5 


(4.04-10 4 ) 


87.5 


Ol(p) 


-0.36 


(-0.32) 


-0.58 


(-0.58) 


-0.19 



The n — p ratio, 

Rnp = ^n/^pi (15) 

the ratio of the parity violating (PV) v.s. the parity conserving (PC) contributions, 

_ bj+ 3(p p + eg 
VP 4 + 3(0* + [Lb) 

Vn = (17) 

and the decay asymmetry parameter, ai, are also given in Table 4. We find that the n — p 
ratio for the direct quark process is much larger than that obtained in the meson-exchange 
calculation. The pure AI = | calculation also yields large R np , indicating that the R np 
enhancement is not related to the AI = | rule violation. 

The results of the one-pion exchange transition are also shown in Tables 3 and 4, for 
comparison |H]. These amplitudes satisfy a n /a p = b n /b p = /„/ f p = y/2, because AI = | is 
assumed for the weak pion-baryon vertex. One sees that the amplitude d in this case is 
dominant. This comes from the tensor part of the one pion exchange and is enhanced due 
to a large relative momentum in the final state. Because the amplitude d is not allowed 
for the nA — > nn by the Pauli principle, the n — p ratio in the pion exchange amplitudes 
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becomes very small. Remember that the tensor part of the direct quark interaction is 
neglected because it is of order {p/m) 2 . 

The magnitudes of the transition amplitudes are in general larger for the meson ex- 
change mechanism than the direct-quark process. It is noticed, however, that some am- 
plitudes, such as a and c, have comparable direct-quark amplitudes. Therefore, if one can 
select the initial and/or final spin states in the decay experiments, it will be possible to 
detect the contribution of the direct quark processes. 

The decay asymmetry parameter describes the angular distribution of the outgoing 
two nucleons in the rest frame, 

W(6) = 1 + a 1 V A P 1 (cos6) (18) 

where V\ is the polarization of A in the nucleus. The parameter a\ is given in terms of 
the two-body decay amplitudes by 

= 2V3(V2c + d)f 

1 a 2 + b 2 + 3(c 2 + d 2 + e 2 + f 2 ) { } 

and thus indicates the interference between the PV and PC components. For the An — > 

nn decay, c and d vanish due to the isospin conservation and therefore a\(n) is zero. 

Recent experiment done at KEK indicates a large negative a\(p) for light hyp er nuclei fl2 



The data is consistent with ai(p) < —0.6. Our calculation yields the correct sign, but 
the magnitude is smaller. The meson exchange calculation done by Ramos et ai, also 
predicted a small magnitude with the correct sign]!]]. 

Recently, Schumacher suggested that analyses of the non-mesonic decays of the 
A=4 and 5 hypernuclei may be usuful in checking the AI = | rule for the nonmesonic 
weak decay. One can parametrize the n-p ratios of the decay of ^He and ^He and the 
ratio of non-mesonic decay widths of ^He and in terms of R n s and R p s, where Rns 
stands for the decay rate of AA" with the initial spin 5|TTJ. Then, using experimental 
data, the ratios of Rns can be extracted. In fact, the R n o/R p o ratio will be 2 for the pure 
AI = | transition, while it becomes 1/2 for the pure AI = 3/2 transition. Because these 
hypernuclei involve only the relative s-wave states, the direct quark processes can play 
significant roles. It may then be possible to see the AI = 3/2 decay, if any. Table 5 shows 
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Table 5: Observables for light hypernuclei. The numbers without parenthesis (in paren- 
thesis) are the results for the parameter I (II). 





full 


AI = 


| omitted 


OPE 


-RnpdHe) 


0.36 


(0.39) 





(0) 


0.09 




0.99 


(0.89) 


1.41 


(1.40) 


0.13 


r n .m. (l^e) /r„ m (^h) 


0.78 


(0.78) 


0.71 


(0.72) 


6.26 


RnO/ RpO 


0.48 


(0.44) 


2 


(2) 


2 



the J = (S — 0) part of the n-p ratio, R n o/R p0 as well as the above observables for our 
direct quark amplitudes with and without the AI = | components. Our full amplitudes 
indeed predict a small ratio R n o/R p o, which indicated a large AI = 3/2 contribution. It is 
therefore concluded that the AI = | components of the direct quark processes is possibly 
observed in non-mesonic decays of light hypernuclei. 

4. Conclusion and Discussion 

We present a quark model calculation of the direct-quark processes of the weak AiV — > 
NN decay, which can be observed exclusively in decays of hypernuclei. We find that the 
transition amplitudes in some of the decay channels are comparable to those in the meson 
exchange decays in magnitudes and show qualitatively distinctive properties. This is 
encouraging because the direct-quark processes may resolve the discrepancies between 
experiment and theory based on the meson-exchange mechanism. Indeed, after averaging 
over the initial and final spin states, we obtain a large n — p ratio from the direct-quark 
amplitudes. In the actual hypernuclear decays, one expects a variety of the spin-isospin 
combinations. We also have to note that the initial orbital angular momentum greater 
than zero may contribute significantly. We, therefore, cannot make a definite prediction 
here. Further study of the AiV — > AW decay with higher partial waves and realistic 
two-baryon wave functions is under way. 
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It is also important to combine the meson-exchange amplitudes with the direct-quark 
ones to make a final quantitative conclusion. We here do not superpose the direct-quark 
and the one-pion-exchange amplitudes because both of them are yet incomplete. The 
one-pion exchange is far from realistic because other mesons, such as K and p, are known 
to contribute to this process significantly. It is not clear whether the phenomenological 
lagrangian for the mesonic decay of A is consistent with the weak effective hamiltonian 
used in the direct-quark calculation. It is important to understand the mesonic decay and 
its A J = \ rule starting from the effective four-quark Lagrangian A study along this 
line is being carried out. 

Our plan also includes the study of AA" — > AW with more realistic wave functions for 
the initial AA" and the final AW states. Especially, effects of the baryon-baryon short- 
range correlation may change the results quantitatively. A more realistic calculation is 
under way. 
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